Casimir forces between arbitrary compact objects: 
Scalar and electromagnetic field 



T. and R. L. JafFe=^ 

^ Institut fiir Theoretische Physik, Universitat zu Koln, Ziilpicher Strasse 77, 

50937 Koln, Germany 

^ Laboratoire de Physique Theorique et Modeles Statistiques, CNRS UMR 8626, 
Bat. 100, Universite Paris-Sud, 91405 Orsay cedex. Prance 

Center for Theoretical Physics, Laboratory for Nuclear Science, and Department of 
Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 

Abstract. We develop an exact method for computing the Casimir energy between 
arbitrary compact objects, both with boundary conditions for a scalar field and 
dielectrics or perfect conductors for the electromagnetic field. The energy is obtained 
as an interaction between multipoles, generated by quantum source or current 
fluctuations. The objects' shape and composition enter only through their scattering 
matrices. The result is exact when all multipoles are included, and converges rapidly. 
A low frequency expansion yields the energy as a series in the ratio of the objects' 
size to their separation. As examples, we obtain this series for two spheres with Robin 
boundary conditions for a scalar field and dielectric spheres for the electromagnetic 
field. The full interaction at all separations is obtained for spheres with Robin boundary 
conditions and for perfectly conducting spheres. 
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1. Introduction 

Casimir forces arise when the quantum fluctuations of a scalar, vector, or even fermion 
field are modified by the presence of static or slowly changing external objects [1]. The 
objects can be modeled by boundary conditions that they place on the fluctuating field 
(f), by an external field, a, to which couples [2], or, in the case of electromagnetism, 
by a material with space and frequency dependent dielectric and magnetic properties. 
The Casimir energy is the difference between the energy of the fluctuating field when 
then objects arc present and when the objects are removed to infinite separation. 

The advent of precision experimental measurements of Casimir forces [3] and the 
possibility that they can be applied to nanoscale electromechanical devices [4, 5] has 
stimulated interest in developing a practical way to calculate the dependence of Casimir 
energies on the shapes of the objects. Many geometries have been analyzed over the 
years, but the case of compact objects has proved rather difficult. In a recent Letter 



Casimir forces between arbitrary compact objects: Scalar and electromagnetic field 2 

we described a new method that makes possible accurate and efficient calculations of 
Casimir forces and torques between any number of compact objects |. The method 
applies to electromagnetic fields and dielectrics as well as perfect conductors. It also 
applies to other fields, such as scalar and Dirac, and to any boundary conditions. In 
this approach, the Casimir energy is given in terms of the fluctuating fleld's scattering 
amplitudes from the individual objects, which encode the effects of the shape and 
boundary conditions. The scattering amplitudes are known analytically in some cases 
and numerically in others. If the scattering amplitudes are known, then the method 
can be applied from asymptotically large separation down to separations that are 
a small fraction of the dimension of the objects. Results at large separations are 
obtained using low frequency and low angular momentum expansions of scattering 
amplitudes. The coefficients multiplying the successive orders in inverse separation 
can be identified with increasingly detailed characteristics of the objects. At small 
separations the manipulation of large matrices, whose dimensions grow with angular 
momentum, eventually slows down the calculation. However at these distances other 
methods, notably the "proximity force approximation" (PEA), apply. Thus it is now 
possible to obtain an understanding of Casimir forces and torques at all separations for 
compact objects. 

The aim of this talk is to provide a pedagogical introduction to our methods 
by treating in detail the simplest scalar field obeying a boundary condition 

on a sharp surface [7]. The complications of electromagnctism and smoothly varying 
dielectrics were already introduced in Ref. [6] and arc discussed here briefiy. Our result 
for the Casimir energy is remarkably simple. For a complex scalar field in the presence 
of two objects it takes the form. 



where the determinant is over the partial wave indices on transition matrices and 
translation matrices U"^ and the integral is over k, = —iuj/c, the imaginary wavenumber. 
For a (real) electromagnetic field the result has the same form with an additional factor 
of 1/2 and the appropriate matrices for scattering and propagation of electromagnetic 
waves. In Sections 4 and 5 the usefulness of this result is demonstrated through several 
specific applications for scalar and electromagnetic fields, respectively. 

The force between atoms at asymptotically large distances was computed by 
Casimir and Polder [8] and related to the atoms' polarizabilities. For compact objects, 
such as two spheres, Feinberg and Sucher [9] generalized this work to include magnetic 
effects. Earlier studies of the Casimir force between compact objects include a multiple 
refiection formalism [10], which in principle could be applied to perfect conductors of 
arbitrary shape. A formulation of the Casimir energy of compact objects in terms of their 
scattering matrices, for a scalar field coupled to a dielectric background, is introduced 




(1.1) 



I This presentation is based on work performed in collaboration with N. Graham and M. Kardar. For 
a complete exposition, see Refs. [6] and [7]. 
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in Ref. [11], where it is suggested that it can also be extended to the EM case. 

Recently Gies et al. [12] used numerical methods to evaluate the Casimir force 
between two Dirichlet spheres for a scalar field, over a range of subasymptotic 
separations, and in other open geometries such as a plate and a cylinder [13] or 
finite plates with edges [14]. Bulgac and collaborators [15] applied scattering theory 
methods to the same scalar Dirichlet problem and obtained results over a wide range of 
separations. The only explicit calculations for subasymptotic distances up to now have 
been for a scalar field obeying Dirichlet boundary conditions on two spheres, a sphere 
and a plate [15] and for electromagnetic fields for a plate and a cylinder [16] and two 
perfectly conducting spheres [6]. 



2. Casimir energy from functional integral: Scalcir field 

In this Section we review formalism essential for our work that is based on the functional 
integral approach of Refs. [17, 18, 19]. 

2.1. Functional integral formulation 

We consider a complex quantum field, 0(x, t), which is defined over all space and 
constrained by boundary conditions C on a set of fixed surfaces Eq, for a = 1, 2, . . . , A^, 
but is otherwise non-interacting. We assume that the surfaces are closed and compact 
and refer to their interiors as "objects." Our starting point is the functional integral 
representation for the trace of the propagator, Tr e~'^HcT/h |-2Qj ^ 

rj^^-iHcT/h _ j ^p^i^ giSsM = Z [C] , (2.1) 

where the subscript C denotes the constraints imposed by the boundary conditions. § 
The integral is over all field configurations that obey the boundary conditions and are 
periodic in a time interval T. 5'[0] is the action for a free complex field, 

SE[(t>\ - jy^j - l^*^!') ' (2-2) 

where the x-integration covers all space. || 

The ground state energy can be projected out of the trace in Eq. (2.1) by setting 
T = —iA/c taking the limit A — > oo, 

£o[C] = - lim ^ In (Tr e"^'^^/'^) = - lim ^ In Z[C] , (2.3) 

§ We have used an abbreviated notation for the functional integral. Since ^ is complex / Vcf) should 
be understood as / 'D(j)'D(j)*, and similarly in subsequent functional integrals. 

II Note that ^ is defined and can fluctuate inside the objects bounded by the surfaces Eq,. In this feature 
our formalism departs from some treatments where the field is defined to be strictly zero (for Dirichlet 
boundary conditions) inside the objects. The fluctuations interior to the objects do not depend on the 
separations between them and therefore do not aifect Casimir forces or torques. 
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and the Casimir energy is obtained by subtracting the ground state energy when the 
objects have been removed to infinite separation, 



fir 

£:[C] = - hm '^\n{Z[C]/Z^) . 

A— »oo 7V 



(2.4) 



In the standard formulation, the constraints are implemented by boundary conditions 
on the field at the surfaces The usual choices are Dirichlet, </> = 0, Neumann, 

(9„0 = 0, or mixed (Robin), 0— A9„0 = 0, where dn is the normal derivative pointing 
out of the objects. To be specific, we first consider Dirichlet boundary conditions. The 
extension to the Neumann case is presented in Ref [7]. As noted in the Introduction, the 
only effect of the choice of boundary conditions is to determine which T- matrix appears 
in the functional determinant, Eq. (1.1). 

Since the constraints on are time independent, the integral over 0(x, t) may be 
written as an infinite product of integrals over Fourier components. 



~ oo 

J me = n ['^M^)]c , 



where 



(2.5) 



(2.6) 



0(x,i)= <^n(x)e™/^, 

n=— oo 

and the logarithm of Z becomes a sum, 

lnZ[C]= ln|y"[P0„(x)],exp J (^(^^J ' l"^ M^) 

As T ^ oo, 'Ylin •^^'^ replaced by |^ rffc, where k = 27m/{cT) and (f)n{^) is 
replaced by 0(x, k). Combining the positive and negative /c-integrals gives 

\nZ[C] dk\n!^J [P0(x, A;)]c exp J c^x (A;2|(/.(x, A;)^ - |V0(x, A;)^) 



TT 

cT 

TT 



dk]n3c{k), 



(2.8) 



where 



5c (A:) = / [P0(x,A;)]cexp 



T 



cix(A;Xx, A;)|'-|V0(x, A;)|') 



(2.9) 



is the functional integral at fixed k. 

To extract the Casimir energy, we use T = —iA/c and Wick rotate the /c-integration 
{k = in with K > 0). % Using Eq. (2.4), we obtain. 



he , ^ 

£[C] = / dKln 

TT Jo 



3oo(«/«) 



(2.10) 



^ A more careful treatment of the rotation of the integration contour to the imaginary axis is necessary 
in the presence of bound states. 
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Here 3c (^'«) is given by the Euclidean functional integral, 



T 

h 



dx («;2|0(x,m)|2 + |V0(x,m)|2) 



(2.11) 



It remains to incorporate the constraints directly into the functional integral using 
the methods of Refs. [21, 17]. Working in Minkowski space, we consider the fixed 
frequency functional integral, 3c (^0 (and suppress the label k on the field 0). Following 
Ref. [21, 17], we implement the constraints in the functional integral by means of a 
functional (5-function. For Dirichlet boundary conditions the constraint reads. 



T 



^^x(^;(x)0(x) +C.C.) 



, (2.12) 



j [P0(x)], = j [P0(x)] n / [T^Q.i^)] exp 

a=l 

where the functional integration over is no longer constrained. Other boundary 
conditions can be implemented similarly. In the resulting functional integral, 

N r- / 

Mk)= n / [^^"W] / [^'^Wlexp cix(A;^|0(x)|2-|V0(x)|2) 

a=l ^ 

+ [ ^^x(^*(x)0(x)+c.c.) 



a 
N 



= n / [^^^?a(x)] / [1^0(x)]exp 



T~ 



(2.13) 



the fields fiuctuatc without constraint throughout space and the sources {g^} fluctuate 
on the surfaces. We denote the new "effective action" including both the fields and 
sources by S[(j),g\. 



2.2. Performing the integral over 

We start with the expression for the fixed-/c functional integral, Eq. (2.13). For any 
fixed sources, {fPa}, there is a unique classical field, 0ci[f?], that is the solution to 
5S'[0, f)]/50(x) = 0. The classical theory defined by S'[0, g]., describes a complex scalar 
field coupled to a set of sources on the surfaces, and is a generalization of electrostatics. 
By analogy with electrostatics, the field is continuous throughout space, but its normal 
derivative jumps by ^»q(x) across E^. Indeed, the classical equations of motion that 
follow from 6S/6(f) — are 

(V' + e) 0ci(x) = 0, for X ^ E„, 
A0ci(x) = 0, for X G Eq,, 

A9„0c4 = ^«(x), for X G E«, (2.14) 

where A0 = 0in — 0out and A(9„0 = 9n0|in — 9„0out- The subscripts "in" and "out" refer 
to the field inside and outside the bounding surface E^. As before, all normals point 
out of the compact surfaces. The solution to Eq. (2.14) is unique up to solutions of 
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the homogeneous equations, which we exclude by demanding that (f>c\ vanish when the 
{Qa} = 0. Continuing the analogy with electrostatics, we can write the classical field in 
terms of the free Green's function and the sources, 

0ci(x) = ^ / rfx'^o(x,x,/c)^^(x') , (2.15) 

where the free Green's function is given by 

^?o(x,x',^) = Ji'^Yx^l = ^kJ2ji{kr^)h'i'\kr^)Yim{m*m{^') 

= ikJ2jiikr<)h'i'\kr^)Yim{^')Y:ji) , (2.16) 

Im 

where the notations refer to whichever of r,r' is the smaller (larger). 

To compute the functional integral over (f), we first decompose into the classical 
part given by Eq. (2.15) and a fiuctuating part, 

(/)(x) = (/)el(x)+(5(/)(x). (2.17) 

Then, because the effective action, S, is quadratic in (p, the 6(j) dependent terms are 
independent of (pd, 

Mk) = n/ [^^i?a(x)]e'^^^'[^l j [P50(x)]exp j dx (fc2|<50(x)|2 - |V50(x)p) 

(2.18) 

The classical action can be simplified by using the equations of motion, Eq. (2.14), which 
make it possible to express the action entirely in terms of integrals over the surfaces 

= \Y.I (^a(x)0cl(x) + C.C. ) , (2.19) 

where 0ci(x) is understood to be a functional of the sources Qa- 

The functional integral over 8(f) is independent of the classical field 0ci and defines the 
energy of the unconstrained vacuum fiuctuations of 0. This term is divergent, or, more 
precisely, depends on some unspecified ultraviolet cutoff. However it can be discarded 

because it is independent of the sources and therefore common to 3c and 3oo- Note that 
this result is an explicit demonstration of the contention of Ref. [22] : the Casimir force 
has nothing to do with the vacuum fluctuations of 0, but is instead a consequence of 
the interaction between fluctuating sources in the materials. It is therefore not directly 
relevant to the fluctuations that are conjectured to be associated with the dark energy. 

Prom Eq. (2.15) it is clear that the solution to Eq. (2.14) obeys the superposition 
principle: 0ci(x) is a sum of contributions from each of the sources, 

0ei(x) = 5^0^(x), (2.20) 

13 
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where satisfies Eq. (2.14) with all sources set equal to zero except for g^. So the 
action can be expressed as a double sum over surfaces and over contributions to 0ci 
generated by different objects. This leaves a partition function, 3c (^), of the form 



(2.21) 



to be evaluated. 



2.3. Evaluation of the Classical Action 

The classical action in Eq. (2.21) contains two quahtatively different terms, the 
interaction between different sources, a ^ f3, and the self-interaction of the source 
Qa- Both can be expressed as functions of the multipole moments of the sources on the 
surfaces. 



2.3.1. Interaction terms: a ^ j3 Consider the contribution to the action from the field, 
(pp, generated by the source, g^, integrated over the surface 

Spa^ \ 1 dKa {gl{K^)M^a) + cc. ) , (2.22) 

J Set 

where the subscript a on indicates that the integration runs over coordinates 
measured relative to the origin of object a. The field 0/j(x^), measured relative to 
the origin of object j3, can be represented as an integral over its sources on the surface 
as in Eq. (2.15). Since every point on is outside of a sphere enclosing E^, the 
partial wave representation of Qo simplifies. The coordinate x> is always associated 
with x^ and x< is identified with x|Jj, so Eq. (2.15) can be written 

M^p) = ikJ2 h^\kr(,)YU±p) [ d^pji{kr'p)Y:^{i^p)gp{^p) . (2.23) 

Note that the arguments of the Bessel functions and spherical harmonics are all defined 
relative to the origin Op. In particular, r'^ and x^ are the radial and angular coordinates 
relative to Op corresponding to a point x' on the surface Yup. The integrals over Yup define 
the multipole moments of the source gp, which will be our final quantum variables, 

Qf3,im = / dxpji{krp)Yi*^{kp)gp{xp), (2.24) 

so that 

(ppi^p) = ik^Qp,imhf\krp)Yim{^p). (2.25) 

Im 

The field (f)p viewed from the surface Eq, is a superposition of solutions to the 
Helmholtz equation that are regular at the origin Oa- Using translation formulas. 
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summarized in Ref. [23], the field generated by object E^j can be written as function of 
the coordinate Xq, measured from the origin Oa, as 

0/3(Xa) = ikY Qp,lmY^l^^'lm3l'{k^a)yVm'{^o) ■ (2-26) 
Im I'm' 

The matrices are shape and boundary condition independent and represent the 

interaction between the multipoles. This result, in turn, can be substituted into the 
contribution Sj^a to the action, leading to the simple result 

SpaiQa, = f E Q*a,l'm'^?^lmQpM + C-C ■ (2.27) 

Iml'm' 

Note that the contributions to the action that couple fields and sources on different 
objects make no reference to the particular boundary conditions that characterize the 
Casimir problem. They depend only on the multipole moments of the fields and on the 
geometry through the translation matrix U"^. 

2.3.2. Self-interaction terms We turn to the terms in where the field and the source 
both refer to the same surface, Sq: 

Sa[Qa] = 1 [ (Q*a(^)M^) + C-C ) ■ (2.28) 

For the self-interactions terms, we only use the coordinate system with origin Oa inside 
the surface 11^, and hence drop the label a on the coordinates in this section. Since 
0a(x) is continuous across the surface, we can regard the (j)^ in Eq. (2.28) as the field 
inside Sq, 0in,a, which is a solution to Helmholtz's equation that must be regular at the 
origin O^, 

0m,a(x) = Y 4)a,lm3l{kr)Ylm{^) ■ (2.29) 
Im 

Substituting this expansion into Eq. (2.28), we obtain 

SalQa] ^IYI {^-mQIm + c.c. ) , (2.30) 

Im 

where the Qa,im are the multipole moments of the sources, defined in the previous 
subsection. 

Finally we relate (j)a.im back to the multipole moments of the source to get an action 
entirely in terms of the Qa,im- The field 0a,out at points outside of obeys Helmholtz's 
equation and must equal (f)a,m on the surface S. Therefore it can be written as (paM plus 
a superposition of the regular solutions to the Helmholtz equation that vanish on Eq, 

0a,out(x) = 0a,m(x) + A0a(x) (2.31) 

= 0a,in(x) +^Xa,/m bK^O>^m(x) + ^ V^',^(A;)4'^(A:r)Y,,„,(x) 
Im \ I'm' 
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The second term, A^^, vanishes on Eq, because T° is the scattering amphtude for the 
Dirichlet problem. 

The field we seek is generated in response to the sources and therefore falls 
exponentially (for k with positive imaginary part) as r ^ oo. Therefore the terms 
in Eq. (2.31) that are proportional to ji{kr) must cancel. Comparing Eq. (2.31) with 
Eq. (2.29), we conclude that Xa,im = —(l>a,imj therefore 



0a,out(x) ^ -J2<f>c.MT.'^lKn'lrnmP{kr)Yi,^,{K) . (2.32) 
Im I'm' 



On the other hand, 0Q,out(x) can be expressed as an integral over the source as in 
Eq. (2.15), 

0a,out(x)= f o?x'^?o(x,x',A;)^«(x'). (2.33) 
Using the partial wave expansion for the free Green's function, Eq. (2.16), we find 

0a,out 

(x) = ikJ2Qc.,i'm'h\}\kr)Y,m'(^), (2.34) 

I'm' 

and comparing with Eq. (2.32), we see that 

ikQa,l'm' — — %,"m'lmi^)^a,lm, 
Im 

or 

0a,/m = -ik Xl[^"]/mVm'<3a,/'m', (2-35) 
I'm' 

where [T"]~^ is the inverse of the Dirichlet transition matrix T"^. When this is combined 
with Eq. (2.30), we obtain the desired expression for the self- interaction contribution to 
the action, 

Sa[Qa] = - Y E Qhm['^''\7ml'm'Qa,l'm' + CC. . (2.36) 
Iml'm' 

2.4- Evaluation of the Integral over Sources 

Combining Eq. (2.36) with Eq. (2.27), we obtain an expression for the action that is 
a quadratic functional of the multipolc moments of the sources on the surfaces. The 
functional integral Eq. (2.21) can be evaluated by changing variables from the sources, 
{Qa} to the multipole moments. The functional determinant that results from this 
change of variables can be discarded because it is a common factor which cancels between 
3c and 3oo- To compute the functional integral we analytically continue to imaginary 
frequency, k — in, k > 0, 

N ( 
a=l'' y a 

(2.37) 
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where we have suppressed the partial wave indices. The functional integral Eq. (2.37) 
yields the inverse determinant of a matrix M"'^ that is composed of the inverse transition 
matrices [T"]~^ on its diagonal and the translation matrices U"^ on the off-diagonals: 

Mf = [TT'Sa3 - U"^(l - S^3) . (2.38) 
Finally we substitute into Eq. (2.10) to obtain the Casimir energy, 

^r^i , ^ detMc(m) 

S[C]^— dKln , ^„ \ { , 2.39 
^ ^ TT Jo detMoo(i«) ' ^ ^ 

where the determinant is taken with respect to the partial wave indices and the object 
indices a, /3, and = [T"]~^5q^ is the result of removing the objects to infinite 

separation, where the interaction effects vanish. 

In the special case of two interacting objects Eq. (2.39) simplifies to 

he r°° 

E^[C] = — / d«;lndet(l - T^U^^W^) , (2.40) 
Ti" Jo 

where T",q; = 1,2, and U"^ are the transition and translation matrices for the two 
objects. 



3. Casimir energy from functional integral: Electromagnetic field 

In this section we provide a brief description of the generalization of the concepts 
from the previous section to electromagnetic fields [6]. Here we start directly with 
a formulation in terms of sources, the current and charge densities J, g. Using that the 
EM gauge and scalar potential [^(x, t), <l>(x, t)] are given by 

[A(x),$(x)] = J dx'Go(x,x')[J(x'),i?(x')] , (3.1) 

where now G'o(x,x') = e^^'l''-^'l/(47r|x - x'|), the action S[J] = J {dk/A7r){Sk[3] + ^^[J]) 
for the currents densities, defined inside the objects, can be written as 

where ^o(x, x') = Go(x, x') — pV ® V'G'o(x, x') is the tensor Green's function. This is 
the analogous expression to the one for a scalar field in Eq. (2.21) with the solution of 
Eq. (2.15) substituted. Next we must constrain the currents to be induced sources that 
depend on shape and material of the objects. Formally this is achieved by integrating 
over currents, inserting constraints to ensure that the currents in vacuum simulate 
the correct induction of microscopic polarization and magnetization (from all 
multipoles) inside the dielectric objects in response to an incident wave. 



Casimir forces between arbitrary compact objects: Scalar and electromagnetic field 11 



Let us consider one object. First, the induced current is Jq = —ikT*a + V x Mq,, 
and since Pq = (ec — 1)E, M„ = (1 — 1//Xq,)B, it can be expressed in terms of the total 
fields E, B inside the object as 

Ja = -ikiea - l)E + V X [(I - l///a)B] . (3.3) 

Second, the total field inside the object must consist of the field generated by and 
the incident field Eo({Ja, x) that has to impinge on the object to induce Jq,, so 
that 

E(x) = Eo({Ja, §"}, ^)+ikJ dx' 6;o(x, x') J,(x') . (3.4) 

The incident field depends on the current density to be induced and on the scattering 
matrix of the object, which connects the incident wave to the scattered wave. It 
is fully specified by the multipole moments of Jq, (see below for details). Substituting 
Eq. (3.4) and B = {l/ik)'V x E into Eq. (3.3) yields a self-consistency condition that 
constrains the current 3 a- If one writes this condition as Cq,[Jq,] = for each object, the 
functional integration over the currents constrained this way for all objects yields the 
partition function 

Z = / n n <^(Ca[Ja Wl) exp {iS[{Ja}]) . (3.5) 

It is instructive to look at two compact objects at a distance d, measured between 
the (arbitrary) origins Oa inside the objects. In this case the action of Eq. (3.2) is 



Sk[{Ja}] = \Y.j t^Xa J;(X„)^E^(X, - d«z) (3.6) 
+ ^ XI Jd^adx'^ Ja(Xa) Goi^a, X^) Ja(x^) , 



where we have substituted the electric field Ea{'^a) = '^k J lix^ ^o(xa, x^) Ja(x^) and 
the fields are measured now in local coordinates so that x = Cq, + Xq,, and da — d {—d) 
for a — 1(2). The off-diagonal terms in Eq. (3.6) represent the interaction between 
the currents on the two materials. A natural way to decompose the interaction between 
charges is to use the multipole expansion. For each body we define magnetic and electric 
multipoles as 

QZ.lrn =ljd^a Ja(Xa) V X [^J i{kr a)Y;j±a)] (3.7) 
Qhm^\jd^aJa{^a)VxV X [^a3l{kra)Y{'^{^a)\ , 

for / > 1, \m\ < /, where A = y^/(r+T), ji are spherical Bessel functions and 
Yim spherical harmonics. We change variables from currents to multipoles in the 
functional integral and, as the final step in our quantization, integrate over all multipole 
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fluctuations on the two objects weighted by the effective action, 

Iml'm' a=l,2 

(3.8) 

with Qf^ = (Qm Zm; Zm)- Formally, this action resembles that for a scalar field. The 
vector nature of the EM field is reflected by the presence of two different polarizations 
and corresponding electric and magnetic multipoles. Let us discuss the terms appearing 
in Eq. (3.8) and sketch its derivation. 

Off-diagonal terms — We need to know the electric fields in Eq. (3.6) exterior to 
the source that generates them. They can be represented in terms of the multipoles 
as E^(x^) = ~k^i^Qi^^°^{xp) where *°^*(x^) are outgoing vector solutions of the 
Helmholtz equation in the coordinates of object /3 +. We would hke to express the 
currents J* in Eq. (3.6) also in terms of multipoles. The difficulty in doing so is that the 
electric field is expressed in terms of outgoing partial waves in the coordinates of object 
(3, while according to Eq. (3.7), the multipoles involve partial waves ^'[^(xq,) that are 
regular at the origin Oa, in the coordinates of object a. Going from the outgoing to 
the regular vector solutions and changing the coordinate system involves a translation 
and change of basis which can be expressed as *°^'(xq, ± dz) = J2i'm' ^UmHm^Tm'i'^oc) 
where the universal (shape and material independent) translation matrices U^^ and 
U^^ represent the interaction between the multipoles. For fixed (Im), {I'm'), they are 
2x2 matrices (magnetic and electric multipoles), and functions of kd only. Their 
explicit form is known but not provided here to save space [24]; they fall off with kd 
according to classical expectations for the EM field. Then the electric field becomes 
^E^(x„ ± dz) = E/m0L*[m(xa) with 0f„ = « Ezm f^imS'm'QL' > ^ud the integration 
in Eq. (3.6) leads, using Eq. (3.7), to the off-diagonal terms in Eq. (3.8). 

Diagonal terms — The self-action, given by the second term of Eq. (3.6), is more 
interesting and more challenging. It can be expressed in terms of multipoles if we use 
the constraint for the currents, Eqs. (3.3) and (3.4). To do so, we first note that in 
scattering theory one usually knows the incident solution and would like to find the 
outgoing scattered solution. They are related by the S'-matrix. Here the situation is 
slightly different. Wc seek to relate a regular solution Eo(xo,) = ^A; ^o,/m^[m (^") 
and the outgoing scattered solution, Ea(xQ) = —kJ2imQTm^"mi'^a), generated by the 
currents in the material — a relation determined by the T-matrix, T" = (§" — I)/ 2 — 
schematically iQ" = T"0q * [25]. We face the inverse problem of determining 0o,;m for 
known scattering data Qf^, hence, 

<t>OM-iT.i'^'']imi'm'Qt'm' (3-9) 

I'm' 

The two components of ^if^^ are given by 1/A; times the weights for E- and M-multipoles of Eq. (3.7) 
with replaced by Yim- Similarly, *5m* have the same expressions upon substituting Bessel by 

Hankcl functions, ji h[^\ 

* Our relation between S- and T-matrix follows Ref. [25] and hence deviates from usual conventions 
by a factor i. 
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so that the incident field is given in terms of the S-matrix, as indicated in Eq. (3.4). Next, 
we express the self-action of the currents inside a body (the second term of Eq. (3.6)), 
as S^[3a] = IJ rfx„[ED* - BH* - (EqD* - BqH*)], the change of the field action 
that results from placing the body into the fixed (regular) incident field Eq = Dq, 
Hq = Bo, where E, H and D, B are the new total fields and fiuxes in the presence of 
the body. Using D = e^E, H = /x~^B inside the body and Eq. (3.3), straightforward 
manipulations lead to the simple self-action <S'^[Ja] = J c?XQ,J*Eo({Ja, §"}). If we 
substitute the regular wave expansion for Eq with coefficients of Eq. (3.9) and integrate 
by using Eq. (3.7), we get Eq. (3.8). 

Having established the action for electric and magnetic multipoles, Eq. (3.7), the 
Casimir energy follows in complete analogy to the scalar case by integrating over the 
multipoles. Hence, the Casimir energy for two objects is again given by Eq. (2.40) with 
an additional factor of 1/2 since the electromagnetic field is real valued. 

4. Applications: Scalar field 

In this section we give a few typical apphcations of our method for a scalar field. We 
consider a real scalar field fiuctuating in the space between two spheres on which Robin 
boundary conditions, — \adn<t> = 0, are imposed. Because a real field has half the 
oscillation modes of a complex field, the Casimir energy in Eq. (2.40) must be divided 
by 2, giving 



We allow for different Robin parameters Ai_2 at the spheres of radius R. This choice 
allows us to study Dirichlet (A/i? — > 0) and Neumann (A/i? — > oo) boundary conditions 
on separate spheres as special cases. We obtain the Casimir energy as a series in R/d 
for large separations d, and numerically at all separations. A comparison of the two 
approaches allows us to measure the rate of convergence of our results. We find that for 
Robin boundary conditions the sign of the force depends on the ratios Xa/R and on the 
separation d. 

4.1. Interaction of two spheres with Robin boundary conditions: general considerations 

The Robin boundary condition — Aq,(9„0 = allows a continuous interpolation between 
Dirichlet and Neumann boundary conditions. Since the radius of the sphere introduces a 
natural length scale, it is convenient to replace A by a dimensionless variable, (a = -^a/ R- 
For > 0, the modulus of the field is suppressed if the surface is approached from the 
outside, while for < it is enhanced. Hence, for negative bound surface states 
can be expected. All the information about the shape of the object and the boundary 
conditions at its surface is provided by the T-matrix. For spherically symmetric objects 
the T-matrix is diagonal and is completely specified by phase shifts Si{k) that do not 




(4.1) 
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depend on m, 

Timi'mik) = Su'S^m'l {e'''^^"^ - l) ■ (4.2) 

In the discussion of the T-matrix for an individual object we again suppress the label 
a. The phase shifts for Robin boundary conditions are 

coUm-"-§fS^, (4.3) 

where ^ = kR and ji [ni) are spherical Bessel functions of first(second) kind. To 
apply Eq. (4.1), we have to evaluate the matrix elements of the transition matrices 
for imaginary frequencies k = m. Using ji{iz) = ^'a/tt/ (2z)Ii^i/2{z) and h\^\iz) = 
—i~^\/2/{'Kz)Kij^il2{z), we obtain for the T-matrix elements 



^ (l/C + l/2)7,+V2(.)-^/,Vv,(.) 

^ ^ ^ 2(i/c+im+v2(^)-^^;-M/2(^)' ^'^ 

where z = kR. 

For two spherical objects we can assume that the center-to-center distance vector is 
parallel to the z-axis. Then the translation matrices simplify. For imaginary frequencies 
the translation matrix elements become 

fl2| 



21 



1, 



ULJ id) = - (-l)-^-''+V(2/ + l)(2/'^-l) 

'/ /' lA 1 1 v 



X y (±i)'"(2r + i, , 

' ^ Mo 




-m 




where d is the separation distance. 

An analysis of the T-matrix shows that it has poles for — 1 < C < 0. For any C, in 
this interval, there exists a finite number of bound states, which increases as ^ 0. In 
the following, we restrict to > and leave the study of interactions in the presence of 
bound states to a future publication. 

The special case of spheres with Dirichlet boundary conditions has been studied in 
Ref. [15]. For two spheres of equal radius, the matrix J^^,, Aj^M^,™? in the notation of 
Ref. [15] is proportional to our T^U^^T^U^i times Ki+ij2{nR)/ Ki,+ii2{nR). It is easy 
to see that this proportionality factor drops out in the final result for the energy if one 
uses Indet = tr In in Eq. (4.1) and expands the logarithm around unity. Thus we agree 
with the results given in Ref. [15]. 

4..2. Asymptotic expansion for large separation 

In this section wc consider the Casimir interaction between two spheres due to a scalar 
field obeying Robin boundary conditions, allowing for a different parameter Ai 2 on each 
sphere. The Casimir energy can be developed in an asymptotic expansion in R/d using 
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Indet = tr In in Eq. (4.1). Expanding the logarithm in powers of N = T^U^^T^U^-^, 
since the T-matrix has no poles in the region of interest we get 



t POO "° 1 

^--^l ^-E-tr(N-). (4.6) 



We have performed the matrix operations using Mathematica. The scaling of the T- 
matrix at small k shows that the p^^ power of N (corresponding to 2p scatterings) 
becomes important at order d^^'^^^^\ Partial waves of order / start to contribute at 
order rf"*^^"*"^') if the T-matrix is diagonal in /, which is the case for spherically symmetric 
objects. Hence the leading terms with p = 1 and / = yield the exact energy to order 
d~^. In the following we will usually restrict the expansion to p < 3, / < 2, yielding the 
interaction to order d"^. 

The large distance expansion of the Casimir energy can be written as 

j=3 ^ ^ 

where bj is the coefficient of the term ~ d~^ . These coefficients can be computed for 
general Robin boundary conditions [7] . Here we restrict to the limiting cases of Dirichlet 
and Neumann boundary conditions. An interesting property is that some coefficients 
bj go to zero for — oo, which corresponds to Neumann boundary conditions. If 
both Xa go to infinity, the coefficients bj vanish for j = 1, . . . , 6, so that the leading 
term in the Casimir energy is ~ d~'^ with a negative amplitude. Hence, Neumann 
boundary conditions lead to an attractive Casimir-Polder power law, as is known from 
electromagnetic field fluctuations. This result can be understood from the absence of 
low-frequency s-waves for Neumann boundary conditions. It is clearly reflected by the 
low frequency expansion of the T-matrix, which has a vanishing amplitude for — > oo 
if Z = 0. If one A^ = and the other goes to inflnity, only the coefficients 63 and 64 
vanish so that the energy scales as d~^ with a positive amplitude. In general, one has 
63 < for Aq, > 0, so that at asymptotic distances the Casimir force is attractive for 
all non- negative finite A^, and for A„ both infinite. It is repulsive if one A^ is finite and 
the other infinite, i.e., if one sphere obeys Neumann boundary conditions. However, at 
smaller distances the interaction can change sign depending on A^, as shown below. 

More precisely, we obtain the following results. If both A^ = 0, the field obeys 
Dirichlet conditions at the two spheres and the flrst six coefficients are 

1 , 1 , 77 , 25 , 29837 , 6491 

03 = — , 04 = — , 05 = , Ofi = , 07 = , bg — . (4.8) 

^ 4' 4' ^ 48' 16' 2880 ' 1152 ^ ^ 

If Neumann conditions are imposed on both surfaces, the coefficients are 



0, 64 = 0, 65 = 0, be = 0, 67 = — — , ^8 = 0, 69 = — uio 



161 3011 _ 175 

96"' ^^-°' ^^-"192" ^^°-"T28' 

(4.9) 
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clearly showing that the asymptotic interaction has a Casimir-Polder power law 
~ 0{d~'^). Also, as in the electromagnetic case, the next to leading order 0{d~^) 
vanishes [6]. Therefore we have included the two next terms of the series. For mixed 
Dirichlet/Neumann boundary conditions, we obtain 

I. n I. n . 17 ^ 11 ^ 663 ^ 235 ,^ 
bs = 0, b, = 0, b, = -, 6e = -, 6.= -, bs = -. (4.10) 

It is important to note that the series in Eq. (4.7) is an asymptotic series and therefore 
cannot be used to obtain the interaction at short distances. 

4-3. Numerical results for Robin boundary conditions on two spheres at all separations 

The primary application of our analysis is to compute the Casimir energy and force to 
high accuracy over a broad range of distances. However, to obtain the interaction at all 
distances, Eq. (4.1) has to be evaluated numerically. We shall see that the domain where 
our method is least accurate is when the two surfaces approach one another. That is 
the regime where semiclassical methods like the proximity force approximation (PEA) 
become exact. Because of its role in this hmit, and because it is often used (with little 
justification) over wide ranges of separations, it is important to compare our calculations 
with the PEA predictions. 

In the proximity force approximation, the energy is obtained as an integral over 
infinitesimal parallel surface elements at their local distance L, measured perpendicular 
to a surface S that can be one of the two surfaces of the objects, or an auxiliary surface 
placed between the objects. The PEA approximation for the energy is then given by 

SpFA=j I S\\{L)dS, (4.11) 

where £\\{L)/A is the energy per area for two parallel plates with distance L. The 
Casimir energy for parallel plates with Robin boundary conditions has been obtained as 
function of A^/L in Ref. [7]. The behavior of the PEA at asymptotically small or large 
Aq/L determines the Casimir interaction as L — > 0. Eor all non-zero values of Ai^2, we 
take \a/L — >■ cxd, but for the Dirichlct case, A = 0, the limit \a/L applies. Eor 
parallel plates with Robin boundary conditions, in the hmit X\^2/L ^ 1 we obtain the 
result for Neumann boundary conditions on both plates, 

$(Ai/L,A2/L)^$o- = -^, (4.12) 

and for Xx^ijL <S 1 we obtain the identical result for plates with Dirichlet conditions. 
Einally for Ai^2/-^ <S 1 ^ A2,i/iv, we obtain the parallel plate result for unlike 
(Dirichlet/Neumann) boundary conditions, 

$(Ai/L, A2/L) ^ = = ^ . (4.13) 
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The last case is relevant at short distances if one of the Aq, = 0. For two spheres of 
radius R and center-to-center separation d with Robin boundary conditions, the PFA 
results can now be obtained easily from Eq. (4.11). In terms of the surface-to-surface 
distance L = d — 2R, we get 



^PFA 



TT 



Rhc 



2{d- 2RY ' 



(4.14) 



where the -|- applies if one and only one = 0, and the — in all other cases. Hence, 
at small separation the interaction becomes independent of Aq, in the sense that it only 
depends on whether one A„ is zero. 

With the results obtained above, we can analyze the sign of the interaction between 
plates and spheres at both asymptotically large and small distances. Since the PFA 
result is expected to hold in the limit where the distance tends to zero, Eq. (4.14) 
predicts the sign of the interaction between spheres in the limit of vanishing distance. 
In the limit of large distances, we can compare the results for parallel planes from 
Eqs. (4.12) and (4.13) to our calculations for two spheres. We find that the sign of 
the asymptotic interaction depends on the choice for Aq, and is identical for plates and 
spheres. Hence, we obtain a complete characterization of the sign of the interaction at 
asymptotically large and small distances for the plate and sphere geometry, which is 
summarized in Table 1. However, as we have seen above, the power law decay at large 
distance is quite different for plates and spheres. 



Ai 


A2 





L — > oo 


remark 












— for all L 


oo 





+ 


+ 


+ for all L 


oo 


oo 






— for all L 


]0,oo[ 


]0,oo[ 






-\- at intermediate L for 

large enough ratio of Ai, A2. 

(for plates: A1/A2 or A2/A1 > 2.8) 


]0,oo[ 





+ 






]0,oo[ 


oo 




+ 





Table 1. The sign of the Casimir force between two plates and two spheres with Robin 
boundary conditions at asymptotically small and large surface-to-surface distance L. 
The sign in these two limits is identical for plates and spheres. Here "— " and 
indicate attractive and repulsive forces, respectively. 



4.. 3.1. Casimir forces for all separations To go beyond the analytic large distance 
expansion, we compute numerically the interaction between two spheres of the same 
radius R with Robin boundary conditions. Guided by the classification of the Casimir 
force according to its sign at small, intermediate and large separations, we discuss the 
six different cases listed in Table 1. Our numerical approach starts from Eq. (4.1). 
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Using the matrix elements of Eq. (4.4) and Eq. (4.5), we compute the determinant and 
the integral over imaginary frequency numerically. We truncate the matrices at a finite 
multipole order / so that they have dimension (1 + Z)^ x (1 + /)^, yielding a series of 
estimates S'^''^ for the Casimir energy. 

£W gives the exact result for asymptotically large separations, while for decreasing 
separations an increasing number of multipoles has to be included. The exact Casimir 
energy at all separations is obtained by extrapolating the series {S'^^^j to Z — > oo. 
We observe an exponentially fast convergence as — ~ Q-s{d/R-2)i^ where 5 is 
a constant of order unity. Hence, as the surfaces approach each other for d — > 2R, the 
rate of convergence tends to zero. However, we find that the first I — 20 elements of 
the series are sufficient to obtain accurate results for the energy at a separation with 
R/d = 0.48, corresponding to a surface-to-surface distance of the spheres of L = 0.083i?, 
i.e., approximately 4% of the sphere diameter. In principle our approach can be extended 
to even smaller separations by including higher order multipoles. However, at such small 
separations semiclassical approximations like the PFA start to become accurate and can 
be also used. 

The results for Dirichlet and Neumann boundary conditions are shown in Fig. 1. 
All energies are divided by £^pfa, given in Eq. (4.14), with the corresponding amplitude 
$0^ (repulsive at small separations) or $q (attractive at small separations). For like 
boundary conditions, either Dirichlet or Neumann, the interaction is attractive at all 
separations, but for unlike boundary conditions it is repulsive. At large separations 
the numerical results show excellent agreement with the asymptotic expansion derived 
above. Note that the reduction of the energy compared to the PFA estimate at large 
distances depends strongly on the boundary conditions, showing the different power 
laws at asymptotically large separations. In the limit of a vanishing surface-to-surface 
distance {R/d ^ 1/2), the energy approaches the PFA estimate in all cases. Generically, 
the PFA overestimates the energy: £pfa is approached from below for R/d 1/2, 
except in the case of Dirichlet boundary conditions on both spheres, where the PFA 
underestimates the actual energy in a range of 0.3 < R/d < 1/2. The deviations from 
the PFA are most pronounced for Neumann boundary conditions. At a surface-to- 
surface distance of L = 3R {R/d = 0.2), the PFA overestimates the energy by a factor 
of 100. 

Casimir interactions for Robin boundary conditions with finite are shown in 
Fig. 2. If Ai = A2 the interaction is always attractive. If the Aq are not equal and 
their ratio is sufficiently large, the Casimir force changes sign either once or twice. 
This behavior resembles the interaction of two plates with Robin boundary conditions. 
However, the criterion for the existence of sign changes in the force now depends not 
only on A1/A2, but on both quantities \i/R and X2/R separately. Even with A1/A2 
fixed, for smaller Xa/R there can be sign changes in the force, while for larger Xa/R the 
force is attractive at all distances. When the ratio A1/A2 is sufficiently large (or formally 
infinite for Dirichlet or Neumann boundary conditions), we can identify three different 
generic cases where sign changes in the force occur: 
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Figure 1. Casimir energy for two spheres of 
radius R and center-to-center distance d: (a) 
Dirichlet boundary conditions for both spheres, 
(b) Neumann boundary conditions for both 
spheres, (c) Spheres with different boundary 
conditions (one Dirichlet, one Neumann). The 
energy is scaled by the PFA estimate of 
Eq. (4.14). The solid curves are obtained by 
extrapolation to I oo. For the smallest 
separation, the extrapolation uncertainty is 
maximal and indicated by an error bar. The 
dashed curves represent the asymptotic large 
distance expansion given in Eq. (4.7) with 
the coefficients of Eqs. (4.8), (4.9) and (4.10), 
respectively. 



First, we consider Dirichlet boundary conditions (Ai = 0) on one sphere and a 
finite non- vanishing X2/R at the other sphere. Figure 2(a) displays the energy for 
X2/R = 10 as a typical example. At large distances the energy is negative, while it 
is positive at short separations with one sign change in between. The asymptotic 
expansion of Eq. (4.7) yields the exact energy at separations well below the sign 
change [7]. While the expansion predicts qualitatively the correct overall behavior 
of the energy, it does not yield the actual position of the sign change correctly. Of 
course, for the Casimir interaction between compact objects, the sign of the force 
JF = —dS/dd is the physically important quantity, not the energy. The distance at 
which the force vanishes cannot be deduced directly from the slope of the curve for 
S = S/SpFA, since one has 



S'id) 



S'id) + 



d-2R 



(4.15) 
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Figure 2. The Casimir energy for two spheres 
with different Robin boundary conditions for 
finite A^: (a) Dirichlet boundary conditions 
and X/R = 10, (b) Neumann boundary 
conditions and X/R = 10, (c) Ai/i? = 10 
and X2/R = 1. The sohd curves correspond 
to extrapolated results for / — > cxd, and the 
dashed curves represent the asymptotic large 
distance expansion given in Eq. (4.7) with the 
coefficients given in Ref. [7]. For logarithmic 
plotting, the modulus of the energy is shown, 
and the sign of the energy is indicated at 
the bottom. The range of separations with 
a repulsive force is shaded. The points 
of vanishing force occur where an auxiliary 
function (dotted curves) is tangent to the solid 
curve, see text for details. 




The force vanishes at the distance do if S'{do) = 0, so that 

^'K) = ^^^^K) . (4.16) 

Hence the distance at which the force vanishes is determined by the position do 
where the curve of the auxihary function t{d) = T{d/ R — 2)'^ is tangent to the curve 
of S. The two unknown quantities do and r are then determined by the conditions 
S{do) = t{do) and S'{do) = t'{do). This procedure allows us to obtain the distance 
at which the force vanishes easily, without computing derivatives numerically. The 
tangent segment of the curve for t{d) is shown in Fig. 2( dotted line. From this 

construction we find that at a distance (i_=^+ the force changes from attractive to 
repulsive for decreasing separations. The position c?_=^_|_ corresponds to a minimum 
of the energy and decreases with decreasing X2/R, so that in the limit A2/i? — >^ it 
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approaches the case of two spheres with Dirichlet boundary conditions, where the 
force is always attractive. 

• Second, we study Neumann boundary conditions on one sphere and a finite non- 
vanishing \2/R at the other sphere. As an example we choose again \2/R = 10, as 
shown in Fig. 2(b). The energy is positive at large distances and becomes negative 
at small distances. The asymptotic expansion is found to be valid well below the 
separation where the sign of the energy changes [7]. Hence, the expansion describes 
the behavior of the energy qualitatively, but does not predict the precise position 
of the sign change. The sign change of the force can be obtained by the method 
described above. At a position d+=^-, the force changes from repulsive to attractive 
with decreasing separation and the energy is maximal. A decreasing (increasing) 
X2/R shifts rf+=>- to smaller (larger) separations. This result is consistent with an 
entirely repulsive (attractive) force for Neumann- Dirichlet (Neumann-Neumann) 
boundary conditions. 

• The third case is obtained if both Aq, are finite and non-zero. A typical example with 
Ai/i? = 20 and A2/i? = 1 is shown in Fig. 2(c). The energy is negative both at large 
and small separations but turns positive at intermediate distances. The asymptotic 
expansion applies again at sufficiently large separations beyond the position where 
the energy becomes positive. For values of the ratio A1/A2 that are larger than 
an it!-dependent threshold, the force changes sign twice, so that it is repulsive 
between the separations and The energy has a minimum (maximum) 
at ((i+^_). If Xi/R increases and X2/R decreases, the repulsive region grows 
until eventually the force becomes repulsive at all separations, corresponding to the 
limit of Dirichlet/Neumann boundary conditions. Decreasing Xi/R and increasing 
X2/R reduces the interval with repulsion. In this case, first the zeros of the energy 
disappear, leaving negative energy at all distances but still a repulsive region, and 
then the two positions where the force vanishes merge, leaving an entirely attractive 
force. 



5. Applications: Electromagnetic field 



As a specific example for the electromagnetic field, we consider two identical dielectric 
spheres. Due to symmetry, the multipoles are decoupled so that the T-matrix is 
diagonal. 



T, 
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Imlm 



= (-1) 



TT 





Ii+iinz) + 2nzll^,{nz) 


-nli^inz) 


'li^^iz) + 2zl^{z)_ 




Ii-^i{nz) + 2nzl'^_^i {nz) 


-nli^{nz) 


'k,^.{z) + 2zKI^,{z) 
(5.1) 



where the sphere radius is R, z — kR, n = e{iK)n{iK), rj = e{iK)/ iJ,{iK), and Ii+i, 
Ki_^_i are Bessel functions. 7]^^^ is obtained from Eq. (5.1) by interchanging e and 
For all partial waves, the leading low frequency contribution is determined by the static 
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electric multipolc polarizability, af = [(e — l)/(e+(/ + l)//)]i?^^"'"^, and the corresponding 
magnetic polarizability, = [(jj — 1)/{ijL+{1 + 1)/1)]R'^''~^^. Including the next to leading 
terms, the T-matrix has the structure 
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M, 3 



,M ,4 



l(2l + l)\\{2l-l)\\^^^''^ ^^''^ 



+ . 



and T^Z^ is obtained by ^ af, 



7;^. The first terms are '-fil 



-[4 + /x(e// + 



H - 6)]/[5(// + 2)2]i^^ = (4/9) [(a* - 1)/(a* + 2)]^i^^ and 7^3, 7^4 are obtained again 
by the replacement, — > e. Now we can apply our general formula in Eq. (2.40) (with 
a factor of 1/2 and the translation matrices for the electromagnetic field [24]) to two 
dielectric spheres with center-to-center distance d. For simplicity, we restrict to two 
partial waves (/ = 2) and two scatterings {p = 1), which yields the exact Casimir energy 
to order d~^^. Matrix operations are performed with Mathematica, and we find the 
interaction 



£ = 



he 



TT 
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- llaf + 867^2 



m] 



d^ 



(5.2) 



, 5473 + E 

'^[«f(77f4-57f4)+E^M]-l + .. ^ 

where E ^ M indicates terms with exchanged superscripts. The leading term, ~ d~'^, 
has precisely the form of the Casimir- Polder force between two atoms [8], including 
magnetic effects [9]. The higher order terms are new, and provide the first systematic 
result for dielectrics with strong curvature. There is no ~ 1/d^ term. 

The limit of perfect metals follows for e — > 00, — > 0. Then higher orders are easily 
included, yielding an asymptotic series 



ho 

TT 'd7 



E = - — ^ 

n=0 



(5.3) 



where the first 10 coefficients are Cq = 143/16, Ci = 0, C2 = 7947/160, C3 = 2065/32, 
C4 = 27705347/100800, C5 = -55251/64, = 1373212550401/144506880, cj = 
-7583389/320, cg = -2516749144274023/44508119040, cg = 274953589659739/275251200. 
This series is obtained by expanding in powers of M and frequency k,, and does not con- 
verge for any fixed R/d. To obtain the energy at all separations, one has to compute 
Eq. (2.40) without these expansions. This is done as before in the case of a scalar field: 
We truncate the matrix N at a finite multipole order I, and compute the determinant 
and the integral numerically. The result is shown in Fig. 3 for perfect metal spheres. 
Our data indicate again that the energy converges as e-^'^'^f^''^^^ to its exact value at 
/ — s> 00, with 5 ~ 0{1). Our result spans all separations between the Casimir-Polder 
limit for d ^ R, and the proximity force approximation (PFA) for R/d 1/2. At a 
surface-to-surface distance L = 4R/3 {R/d = 0.3), PFA overestimates the energy by a 
factor of 10. Including up to / = 32 and extrapolating based on the exponential fit, we 
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Figure 3. Casimir energy of two metal spheres, divided by the PFA estimate 
fpFA = -{TT^/lUO)hcR/{d - 2Rf, which holds only in the limit R/d 1/2. The 
label I denotes the multipole order of truncation. The curves / — oo are obtained 
by extrapolation. The Casimir-Polder curve is the leading term of Eq. (5.3). Inset: 
Convergence at short separations. 

can accurately determine the Casimir energy down to R/d = 0.49, i.e. L = 0.04/2. A 
similar numerical evaluation can be also applied to dielectrics [7]. 
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